Observation of Distinct Electron-Phonon Couplings in Gated Bilayer Graphene 
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A Raman study of a back gated bilayer graphene sample is presented. The changes in the Fermi 
level induced by charge transfer splits the Raman G-band, hardening its higher component and 
softening the lower one. These two components are associated with the symmetric (S) and anti- 
symmetric vibration (AS) of the atoms in the two layers, the later one becoming Raman active due 
to inversion symmetry breaking. The phonon hardening and softening are explained by considering 
the selective coupling of the S and AS phonons with interband and intraband electron-hole pairs. 

PACS numbers: 63.20.Kd, 78.30.Na, 81.05.Uw 



The interaction of electrons and phonons is a funda- 
mental issue for understanding the physics of graphene, 
resulting in the renormalization of phonon energy and 
the break-down of the adiabatic (Born-Oppenheimer) ap- 
proximation [l|,i,il|. Although many new physical in- 
sights about the electron-phonon interaction were found 
in monolayer graphene 0, B l4, S 16( , the bilayer graphene 
is a unique system to study distinct couplings between 
electrons and phonons that have different symmetries 0] ■ 
In this paper we show that the application of a gate volt- 
age in bilayer graphene splits the symmetric and anti- 
symmetric optical phonon components, confirming a re- 
cent theoretical prediction for the distinct interactions 
of these phonons with intraband and interband electron- 
hole pairs [7|. This result is specially relevant since the 
development of bilayer graphene devices with tunable gap 
depends on the detailed understanding of the interaction 
between electrons and phonons in this material [1, H, [13] ■ 

The zone center E2g phonon mode of monolayer 
graphene, which gives rise to the Raman G-band around 
1580 cm"'^, exhibits a very strong coupling with electron- 
hole pairs. This interaction renormalizes the phonon en- 
ergy hwQ^ giving rise to the so-called Kohn anomaly at the 
zone center of the phonon dispersion 1]. Moreover, theo- 
retical models have predicted an interesting behavior for 
the phonon energy when the Fermi level ei? is changed by 
varying the electron or hole concentration: the phonon 
energy softens logarithmically for values of the chemical 
potential smaller than half of the phonon energy, and it 
hardens otherwise. This hardening is ascribed to the sup- 
pression of the interaction between phonons and electron- 
hole pairs for \/S.£f\ > hujQ/2. Besides this, an increase in 
either electron of hole density increases the phonon life- 
time 0, H, 0, H, Q due to the inhibition of the process of 
phonon decay into electron-hole pairs, thus reducing the 
G-band lincwidth. Recent Raman studies of monolayer 
graphene have confirmed the theoretical prediction con- 
cerning the hardening and narrowing of the the G-band, 
by doping the sample with electrons or holes 0, HI, Q . 

In the case of bilayer graphene, the E2g phonon mode 
splits into two components, associated with the symmet- 
ric (S) and anti-symmetric (AS) displacements of the 
atoms in the two layers. Moreover, due to the splitting 
of the TT and tt* bands in this material, phonons can cou- 
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FIG. 1: (a) Optical microscope image of gated graphene bi- 
layer sample (the black scale bar corresponds to 15 fim). (b) 
Raman spectrum of the G' band of bilayer graphene. (c) Ra- 
man spectrum of the G band for different values of the applied 
gate voltage. 



pie with electron-hole pairs produced in interband or in- 
traband transitions. T. Ando [3] calculated recently the 
self-energy of the S and AS phonons for varying Fermi en- 
ergies and predicted the hardening and softening of the S 
and AS phonons, respectively, induced by electron or hole 
doping. In this work we present experimental results that 
confirm the theoretical prediction by T. Ando fi\, thus 
showing that the G-band of bilayer graphene has indeed 
two components which exhibit opposite dependence as 
the Fermi level is tuned. 

Figure [ifa) shows an optical microscope image of 
the bilayer graphene device obtained by exfoliating bulk 
graphite on a 300 nm thickness silicon oxide layer on 
the top of a heavily p doped silicon substrate. An elec- 
trical contact was made by soldering an indium micro 
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wire directly on the top of the flake 
were induced in the sample by applying a gate voltage 
Vg to the Si substrate with respect to the graphene con- 
tact. The Raman measurements were carried out at room 
temperature using a triple Dilor XY spectrometer with 
a resolution smaller than 1 cm~^, an 80x objective with 
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FIG. 2: (color online) (a) Raman G-band of the bilayer 
graphene for -80 V, -40 V, -20 V and +40 V gate voltages 
. Two Lorentzian curves are needed to fit the G-band for -80 
V, -40 V and -20 V. (b) Displacement of the atoms for the 
symmetric and anti-symmetric highest energy phonon modes 
in the F point of bilayer graphene. 



spot size of ^ 1 /im and the 2.41 eV laser excitation with 
1 mW power. 

The characterization of the number of layers in our 
graphene sample was determined by analyzing the shape 
of the G' band around ~ 2700 cm~^ which is shown in 
Fig. [IJb) (this band is also called as D* or 2D by other 
authors) . The G' band of our sample is composed of four 
peaks and exhibits a typical shape of a bilayer graphene 

[il[il[3[il. 

Figure[T{c) shows the Raman spectra taken at different 
applied gate voltages. We observe in this figure that both 
the position and the shape of the G-band depend on the 
applied gate voltage, in agreement with previous Raman 
studies of gated bilayer graphene [lBi[l3]- However, as 
our experiments were done at room temperature, we were 
not able to observe the initial softening of the G-band for 
Aej?! < fvjjQ/2, as observed at low temperatures by Yan 
et al. 13 • 

Figure [2{a) shows the fittings of the Raman G-band of 
the bilayer graphene device for different values of applied 
gate voltage. The G-bands observed in the spectra taken 
at Vg equals to -80 V, -40 V and -20 V clearly exhibit 
two peaks. All spectra taken at positive values of gate 
voltage could be fitted by a single Lorentzian. 

Figures [S^a) and [3jb) show the relative shift (with re- 
spect to G-band position in the -I- 50 V spectrum) and the 
full width at half maximum (FWHM) of the Lorentzians 
that fit the G-band as a function of the applied gate volt- 
age, respectively. Notice that the frequency of the high 
energy component increases with decreasing values of ap- 
plied gate voltage, whereas the opposite behavior is ob- 
served for the low energy component. Figure [Sjb) shows 
that the FWHM exhibits a maximum for Vg = -1-50 V 
spectrum, and it clearly decreases with decreasing values 
of apphed gate voltages. 

Let us discuss the origin of the two peaks that compose 
the G-band shown in Fig. [Ifa). The bilayer graphene is 
formed by two graphene layers in the AB Bernal stack- 
ing, where the point group at F point is Da^j. The 
phonon branch associated with the E2g mode of mono- 
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FIG. 3: (color online) (a) The black squares and red circles 
correspond to the relative shift (Awq) of the symmetric and 
anti-symmetric components of the G band, respectively, with 
respect to the position in the -1-50 V spectrum, as a function of 
the applied gate voltage (bottom scale) . The solid curves cor- 
responds to the theoretical prediction for the relative G-band 
shift 01 in units of Awo as a function of the Fermi energy (top 
scale), (b) The black squares and red circles correspond to 
the FWHM of the symmetric and anti-symmetric components 
of the G-band, respectively, and the solid curve corresponds 
to the theoretical prediction for the symmetric phonon mode 
0- 



layer graphene gives rise to two branches for bilayer 
graphene, one symmetric (S) and other anti-symmetric 
(AS) components (in-phase and out-of-phase displace- 
ments of the atoms in the two layers), which are rep- 
resented at Fig. [2{b). At the center of the Brillouin zone 
(F point), the symmetric and anti-symmetric vibrations 
belong to the two doubly degenerated representations Eg 
and E„, respectively. The anti-symmetric E„ mode is not 
Raman active since the Y)zd group is centro-symmetric. 
However, if the inversion symmetry operation of the bi- 
layer graphene is broken, the system is described by the 
G-iv point group at the F point. In this case, both the 
symmetric and the anti-symmetric modes belong to the 
E representation, and are Raman active. 

The inversion symmetry breaking can be due to the 
different materials that the top and bottom graphene 
layers of our device are exposed to and/or to a non- 
homogeneous doping of the top and bottom layers. The 
non-equivalence between the top and bottom layers de- 
creases the symmetry, making the anti-symmetric mode 
active in the Raman spectra, explaining thus the two 
components of the G-band. The energy separation be- 
tween the S and AS components is negligible when the 
Fermi level is close to the Dirac point, and this can be 
seen in the spectrum taken at -I- 40 gate voltage in Fig. 
[IJb). However, the doping with electrons or holes induces 
a measurable energy splitting for those components, as 
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shown in the spectra between -20 V to - 80 V in Fig. 

One could argue that the observed spHtting of the G 
band is related to symmetric mode coming from the top 
and bottom layers, in the case of a non-homogeneous dop- 
ing of the two layers. However, this interpretation has 
serious drawbacks in the present case. Notice that the 
low energy component of the G-band is about five times 
less intense than the higher energy one. Considering that 
the transmitance of visible light in monolayer graphene 
is about 97.7% [IS], the contributions of symmetric mode 
from the upper and bottom layers for the Raman scat- 
tering are expected to be almost the same. Furthermore, 
our results show that the energy of the low energy com- 
ponent decreases when the Fermi level is changed. This 
result would not be expected if one associates the lower 
energy Raman peak with the symmetric mode of the bi- 
layer graphene bottom layer. We stress the fact that the 
splitting of the G-band was not reported in the two re- 
cent studies of gated bilayer graphene devices de- 
spite the fact that the G-band observed in reference [lB| 
is clearly asymmetric. Possibly, our observation of this 
splitting is due to the characteristic of our device which 
leads to the inversion symmetry breaking and allows the 
observation of the non-active Raman AS vibration. 

In order to discuss the phonon renormalization effect 
in bilayer graphene, we must consider the selection rules 
for the interaction of the symmetric and anti-symmetric 
phonons with the interbands or intrabands electron-hole 
pairs. The band structure of bilayer graphene near the 
K point shown in Figure 2] consists of four parabolic 
bands, two of them touch each other at the K point, and 
the other two bands separated by 271, where 71 ^ 0.35 
eV. The electron-phonon interaction in bilayer graphene 
is described by a 2x2 matrix for each phonon symme- 
try, where each matrix element gives the contribution 
of electron-hole pairs involving different electronic sub- 
bands For the symmetric phonon mode, all matrix 
elements are different from zero, and this phonon can 
interact with both interband or intraband electron-hole 
pairs, as shown in Figure IJfa), giving rise to the phonon 
energy renormalization (Kohn anomaly). However, for 
the anti-symmetric phonon mode, the diagonal terms of 
the matrix are null, showing that there is no coupling 
between AS phonons and interband electron- hole pairs. 
Therefore, no Kohn anomaly is expected for the antisym- 
metric phonon mode when the Fermi level is at the Dirac 
point. 

However, if the Fermi energy is changed [for instance, 
ep < as shown in Fig. (Hb)], intraband electron- 
hole pairs can be produced by phonons. Now, the anti- 
symmetric phonons have also their energies renormalized, 
giving rise to the Kohn anomaly. Notice that as the en- 
ergy separation between the tti and 7T2 bands (71 0.35 
eV) is larger than the G band energy 0.2 eV), the 
electron-hole pairs creation by those phonons is a vir- 
tual process. Therefore, we expect the dependence of the 
shift of the anti-symmetric mode with Fermi energy to be 



smoother than that of the symmetric mode. A similar re- 
sult has been recently observed in gated semiconducting 
carbon nanotubes [23| which have the band gap larger 
than the phonon energy. 




FIG. 4: (color online) Parabolic band structure of bilayer 
graphene near the K point. The vertical arrows illustrates 
the possible transitions induced by symmetric (green) and 
antisymmetric (red) q=0 phonons for (a) interband electron- 
hole pairs creation at (eF=0) and (b) intraband electron- hole 
pairs creation (ef <0). The gap opening is not considered in 
this diagram. 

In order to explain our experimental results, we have 
also plotted in Figs. Fig. [3] (a) and (b) frequency shift 
and broadening of the G band calculated by T. Ando 0] 
as AujG = i?e[H(±)(w)] and Tq = -Jm[H(±)(w)] respec- 
tively, where H is the phonon self energy given of the 
symmetric (+) or anti-symmetric (-) components, given 
by: 

J 

lfisBjk)~fie^'j>k')]{i^sjk-es'j'k') \ 
(huj+z5y^-(e,,k-e,>,>k'r • ^ ' 

The functions Ssjk are the electronic band dispersions 
{k is the modulus of the wavevector), s = +1 and -1 
denotes the conduction and valence bands, respectively, 
and J = 1, 2 specifies the two bands within the valence 
and conduction bands. The function f{ssik) is the Fermi 
distribution and the matrix elements <l>j^, (fc) gives the 
relative contribution of the S (-I-) and AS (-) phonons for 
the interband and intraband electron-hole pairs formed 
in transitions 0- The damping parameter S describes 
the charge inhomogeneity in the sample and washes out 
the logarithmic singularity present at Ei? = hioo/2. We 
have considered the value of 5 = 0.1 eV in the theoretical 
curves depicted in Fig. [31 

Theoretical curves Awg and Tq are calculated in units 
of Awq, where = 0.196 eV is the energy of the optical 
phonon and A is related to the strength of the electron- 
phonon coupling and given by 

A = 0.16x 10-^(A'el/-2)[^]2 (2) 
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In order to compare the experimental and theoretical 
results depicted in Figs. [21(a) and (b) we need to convert 
the experimental horizontal scale (applied gate voltage, 
in the bottom axis of Fig. [3]) into the theoretical scale 
(Fermi energy, in the upper axis of Fig. [3|), and the 
experimental vertical scales (in cm~^, left side axes of 
Fig. [3]) into the theoretical scales (in units of XllIq, right 
side axes of Fig. 

To match the scales of the experimental (left axis) and 
calculated (right axis) values of Aujq and Fg, shown in 

Fig. [3l we have used the value d^o/db — 6.4 eV A 
in Eq. [5] [l^ . Both theoretical curves Awg and Tq are 
plotted in Fig. [3] as a function of the Fermi energy ep 
(top axis) , while the experimental values are plotted as a 
function of gate voltage (bottom axis). In order to scale 
these axes, we have used the fact that both ep and Vg 
are related to the electron (or hole) density n, which is 
assumed to be the same in both layers. The modulus of 
the Fermi energy [efI is related to the electron (or hole) 
density in the regime lep] < 71 by 0: 



IefI = ^(-71 + V4"^7' + 7?), (3) 

where 7 — ^ajo, a — 2A2A is the lattice constant, 
7o ~ 3eV is the in-plane nearest-neighbor tight-binding 
parameter and 71 ^ 0.35 eV is the out-of-plane nearest- 
neighbor parameter [see Fig. [lib)]. A parallel plate ca- 
pacitor model gives n — 7.2 x 10^°cm~^y~^(Vg - Vd), 
where Vd is the gate voltage needed to move the Fermi 
level to the Dirac point due to intrinsic doping of the 
sample. The value of Vd = -|- 50 V was obtained from 
the spectrum where the G-band has the largest FWHM 
[see Fig. intb)], since the symmetric phonon lifetime is 
minimum for ej? = 

Figure [3] shows that there is a good qualitative agree- 
ment between our experimental results and the theo- 



retical prediction of the dependence of the frequency 
and width of the symmetric and anti-symmetric optical 
phonons in bilayer graphene with the Fermi level. We 
shall stress that this good agreement is obtained despite 
the different approximations that have been done. Notice 
first that the theoretical model doesn't take into account 
the trigonal warping effect, the angular dependence of 
the electron-phonon coupling, and the strength d^i/db 
for the coupling between electrons and AS phonons [see 
Fig. IDJb)]. Moreover, we have used the simple model of 
a parallel capacitor with homogeneous carrier concentra- 
tion in the two layers. A more complete model is needed 
to improve the fitting of the experimental data. 

In summary, we have shown in this work that bilayer 
graphene is a unique material where the phonon renor- 
malization tuned by charge transfer depends strongly on 
the symmetry of the phonons involved in the creation 
of electron-hole pairs. The inversion symmetry breaking 
of our bilayer graphene device allows the observation of 
both the symmetric and anti-symmetric phonon modes 
in the Raman G band. The hardening of the symmet- 
ric mode and the softening of the antisymmetric phonon 
mode when the sample is doped is explained by the se- 
lection rules associated with the creation of interband 
or intrabands electron-hole pairs by symmetric and anti- 
symmetric optical phonons, giving an experimental sup- 
port for the theoretical predictions of the optical phonons 
dependence on charge concentration in bilayer graphene 
0. 
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